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We derive the two-loop corrections to Bhabha scattering from heavy 
fcrmions using dispersion relations. The double-box contributions are ex- 
pressed by three kernel functions. Convoluting the perturbative kernels 
with fermionic threshold functions or with hadronic data allows to deter- 
mine numerical results for small electron mass m e , combined with arbitrary 
values of the fermion mass m/ in the loop, rr? e << s, t, m^, or with hadronic 
insertions. We present numerical results for m/ = m^,m T ,m t op at typical 
small- and large-angle kinematics ranging from 1 GeV to 500 GeV. 

1. Introduction 

Bhabha scattering is one of the theoretically best studied scattering 
processes at e + e~ colliders and can also be measured with a high precision. 
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The accuracy of the Monte-Carlo programs developed originally for physics 
at LEP is about 10 -3 , and with a complete two- loop calculation one may 
reach 10 -4 . The latter number is indicative of efforts for the International 
Linear Collider (ILC), here especially in the GigaZ option running at the 
Z boson resonance, but also for meson factories running at much smaller 
energies of about 1 or 10 GeV. 

Recent years brought considerable progress in the determination of the 
virtual NNLO corrections. The virtual 0(a 2 ) contributions to the massless 
differential Bhabha cross section have been determined in [2]. Shortly after, 
this result was used for deriving the 0{a 2 L) (L = ln(s/m 2 )) corrections to 
massive Bhabha scattering in [3J . The missing photonic correction terms of 
order 0(a 2 L°) were derived, also from [2j, in [HE]. The virtual corrections 
from fermion loop insertions, including the corresponding double-box dia- 
grams, could not be covered by that method. For nj = 1, i.e. the case of 
only electron loops, the corresponding diagrams were evaluated analytically 
in [njEyHyH] , and the net rif = 1 cross section in [TUJ[TT]. 

At this stage the numerically most important two-loop corrections were 
known. For a complete treatment one needs additionally the nj = 2 two- 
loop corrections with heavy fermion insertions, including the hadronic cor- 
rections which replace the loop insertions from light quarks. The leptonic 
rif = 2 contributions have been derived quite recently in two papers in 
the limit m 2 « m 2 « s,t; with a direct Feynman diagram calculation 
in [12], and using a factorization formula that relates massless and massive 
amplitudes in [13] (for that method see also [H]). 

It might be interesting to mention that the original expectations on the 
necessity of a complete, direct two-loop massive Feynman diagram evalua- 
tion were not fulfilled. After the analytical evaluation of a massive planar 
and a massive non-planar double-box diagram (both with seven propaga- 
tors) in [15] and in [16j . resp., there was hope to evaluate all the remaining 
diagrams soon. There are 33 two-loop box master integrals, nine of them 
with seven lines |17j . In fact, from recent studies on the exact and mass ex- 
panded treatment of two-loop box master integrals in [18J and [19], resp., it 
became clear that an evaluation of all the massive master integrals is a more 
complicated task than was expected. Quite recently, the case of non-planar 
master integrals was successfully treated in another, but related context |2U| . 
Proceeding similarly for Bhabha scattering seems feasible now. As a matter 
of fact, due to these reasons, the direct Feynman diagram approach was not 
used for the phenomenologically needed two-loop predictions and the above- 
mentioned papers [15]- [19] remained so far a mere interesting, challenging 
theoretical development. 

In this paper, we report on the evaluation of the leptonic nj = 2 two- 
loop contributions with arbitrary mass of the heavier fermion, i.e. exploring 
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Fig. 1. Classes of two-loop diagrams for Bhabha scattering containing at least one 
fermion loop. 

the extended kinematical region ml « mj,s,t. We use the dispersion 
approach, so that our formulae may be applied without further modification 
also to hadronic corrections. 

For a review on the status of Monte-Carlo studies for Bhabha scattering 
at this conference we refer to [21J, and for a discussion on progress related 
to radiative loop corrections to [2"2"] . 

2. Formulae 

The classes of two- loop nj = 2 corrections are shown in Figure [TJ The 
four direct and four crossed fermionic two-loop box diagrams, obtained ap- 
plying proper permutations to the sample diagram shown in Figure[TJ are in- 
frared (IR) divergent, and they have to be combined with other IR-divergent 
factorizable corrections in order to get an IR-finite contribution to the cross 
section. In particular, we have to add the interference of two-loop box (class 
e) and reducible vertex diagrams (class c) with the tree-level amplitude to 
the interference of one-loop vertex and box diagrams with one-loop vacuum 
polarization diagrams. Finally, we construct an IR-finite quantity taking 
into account also the real emission of one soft photon from one-loop vac- 
uum polarization diagrams. Sample contributions are given in Figure [2j 



The net contribution of pure self-energy corrections (classes a — b), ir- 
reducible vertex diagrams (class d), and the aforementioned IR-divergent 
contributions reads as 

d(J NNLO,ferm. ^ ^a.b d(J d d(J vest 



Fig. 2. The fermionic two- loop boxes combine with other diagrams to an infrared- 
finite cross-section contribution. 



where da 



/dfl can be splitted in two components, 
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(2) 



We concentrate now on the renormalized two-loop box diagrams of class 
e, whose total contribution to the cross section may be written as 




Here the auxiliary functions A s and At can be conveniently expressed through 
three independent form factors B\, with i = A, B, C, evaluated with different 
kinematical arguments, 



The particular contribution of the diagram of Figure Q] coming from the 
interference with the tree-level s-channel is B A (s, t), and from the t-channel 
is B B (s,t). 

For the evaluation of hadronic corrections, we observe that each term 
of Eq. ([!]) can be written through the convolution of a kernel function Kj, 
I = A,B,C, with the hadronic cross-section ratio Rhad, 



A s 
A t 



B A (s,t) + B B (t,s) + B c (u,t) - B B (u,s), 
B B (s,t) + B A (t,s) - B B (u,t) + B c (u,s). 



(4) 




— Rhad{z)K I {s,t,z). 
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Fig. 3. The master integrals for the two-loop box kernel functions. 



For leptons and the top quark, we have to replace 4 M% — > 4 mj and Rhad 
Rfer, given by 




Rfer(z) = Q} Cf \ 1 - 4—^ 1 + 2 + € ^ er (z), (6) 



where is the electric charge in units of |e| and is the color factor. Being 
the box diagrams IR divergent, and showing poles in e, one should take into 
account higher orders in e for R. However, after assembling box diagrams 
with factorizable corrections, IR poles cancel and R can be evaluated at 
order e°. 

The three box kernels are our main technical result. They have been 
derived with the aid of the master integrals of Figure in the limit rr? e < < 
mj, s, t. The master integrals were determined with IdSolver and evaluated 
with the Mathematica packages ambre [23J and MB [24J, and eventually mass 
expanded with a Mathematica package. We also made use of FORM We 
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reproduce one of the kernels here0 



Kc(x,y,z) 



I {2^x 2 L x + AQ 2 x 2 ( Z --2)-2{x 2 + 

v — z) { e Vii / 



3 m 2 . (y — z) 

+ xy)L x + x 2 (~ - l) L y + 2x 2 (J - l) L 2 + Ax 2 L x L y 



+ x 2 ( Z -- l) ln(4) - 2x 2 ( Z - - \) In 2 (4) + Ax 2 (*- 



y i vm^ ^y 2/ ^m^' \y 

1) ln(4) ln(l--)+2x 2 \J—^)L x -x 2 (- + V - 

J Vrrig/ v y' \y z 

2) ln(l-^)-4xMn(l-£)L t + 4^(£-l)Li 2 Q 
2^Lb(l + i)}, (7) 



where L x = ln(— m 2 /x), L y = ln(—m 2 /y) and F e is the normalization factor 



m 2 7r e 7B 



(8) 



Here [i is the 't Hooft mass unit and 7^ is the Euler-Mascheroni constant. 



3. Numerical Results 

The sum of the box contributions with IR-divergent factorizable correc- 
tions (see Figure [2] for sample cases) is infrared-finite and can be cast in the 
following form, 



da 



rest 



1 



IT J S ^ J A M" 2 Z t — Z 



AM 2 

_ f°° , R(z) 1 
+ Re / dz— — 

Jam 2 z s — z + 10 

coo 



Fi(z) 



F 2 (z) + F 3 (z) ln(l - — — 
v s + 10 



R(z) 



+ 7rlm / dz 

Jam 2 z s — z + id 

We may show here one of the auxiliary functions Fi(z), 
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1 Some formulae have to be omitted here due to limited space; they may be found at 
the webpage http://www-zeuthen.desy.de/theory/research/bhabha/bhabha.html. 
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(10) 



Table [U and Table [2] contain numerical results for small- and large-angle 
scattering at a variety of energy scales. We report the QED tree-level pre- 
diction, as well as the process-dependent contributions at NNLO of Eq. ([9]); 
in other words, we exclude from the tables pure self-energy corrections, 
which can be described introducing a running fine-structure constant and 
were deeply investigated in the past (see [26] ) , and irreducible vertex contri- 
butions (see [27] and [28] ). A complete phenomenological analysis requires 
also to add the corresponding terms arising from unresolved real fermion 
pair production. 
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The Standard Model cross sections shown here rely on Born formulae 
with Z boson and photon exchange a la Zfitter [29,30j. Moreover, al- 
though the contributions from electron loops have been obtained by exact 
(in m e ) evaluation of the Feynman diagrams, we report here for consistency 
the approximated results for m e « s,t. 

We further compare the analytical results of [12] with those obtained 
with the dispersion approach and it is nicely seen that the former approach 
the latter in regions where the former are expected to become good ap- 
proximations. In both tables, for each fermion flavour, we show the result 
obtained through the dispersion-based approach (first line) and the one com- 
ing from the analytical expansion (second line), neglecting 0(mj/x), where 
x = s,\t\,\u\. When m 2 > x, the entry is suppressed. We switch off in 
the tables the effect of the logarithm containing the energy of soft photons 
setting uj = y/s/2. 

The heavier fermions have less influence on the net result, and the 
top quark decouples nearly completely. Between 1 and 500 GeV the sum 
of boxes with factorizable diagrams with muon loops contributes, roughly 
speaking, at the order of permille to the net pure QED cross section, and 
the tau lepton contributes less. The Z resonance distortes this figure, by 
making the influence of two-loop contributions less influential for large-angle 
scattering where the resonance dominates. 



[GeV] 


1 


10 


M z 


500 


QED Born 


440994 


4409.94 


53.0348 


1.76398 


rest e 


193 


5.73 


0.1357 


0.00673 


A* 


< 1 


0.42 


0.0408 


0.00288 




X 


0.08 


0.0407 


0.00288 


r 


< 1 


< lO" 2 


0.0027 


0.00088 




X 


X 


-0.0096 


0.00084 


t 


< 1 


< 10" 2 


< 10" 4 


< 10" 5 




X 


X 


X 


X 



Table 1. Numerical values for the differential cross section in nanobarns at a 
scattering angle 8 = 3°, in units of 10 2 ; Mz = 91.1876 GeV. Bold-face entries are 
obtained with dispersion relations. The energy of the soft photon is chosen to be 
uj = s/s/2. 



4. Summary 

We have evaluated the nj = 2 virtual two- loop corrections to Bhabha 
scattering due to fermions with arbitrary mass my in the limit of vanish- 
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lull Born 


466558 


1 /■ill'' • 1 ' — 

4686.27 


1273.2680 


0.85496 


rest e 


807 


14.53 


0.2706 


0.01193 




160 


6.08 


0.1470 


0.00726 




153 


6.08 


0.1470 


0.00726 


T 


2 


1.33 


0.0752 


0.00457 




X 


1.07 


0.0752 


0.00457 


t 


< 1 


< 10" 2 


0.0005 


0.00043 




X 


X 


X 


-0.00013 



Table 2. Numerical values for the differential cross section in nanobarns at a 
scattering angle 9 = 90°, in units of 10 -4 . See Table [T] for further details. 

ing electron mass m e . We have not combined these (infrared finite) virtual 
contributions with those arising from irreducible vertices; the latter are 
logarithmically enhanced by terms of order up to ln 3 (s/m 2 ), but are in- 
dependent of ln(s/mg), being collinear finite. They have to be assembled 
with the unresolved real heavy fermion emission, which is known to cancel 
the ln 3 (s/m 2 ) and might lead to a suppression of the net effect. For the 
phenomenological analysis we will have also to take into account the effect 
of the running of the fine-structure constant. Concerning the results shown 
in the tables, the numerical contributions do not exceed the per mille level, 
and depend strongly on the kinematics. The formulae presented here apply 
to the leptons n and r, but also to the top quark. The latter decouples at 
small energies, but has to be taken into account at the ILC. 

The determination of hadronic corrections is in preparation. 

Note added: 

After completion of this article, a draft [31J appeared, where the authors 
also study the fermionic corrections to Bhabha scattering with arbitrary 
masses of the internal fermions. 
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